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Preface 


This report illustrates the evolution of directionally distributed, shallow- 
water waves at times of high energy as measured with a high-resolution direc- 
tional wave gauge. The work was motivated by the need for a complete de- 
scription of storm wave fields incident on the nearshore to augment under- 
standing and modeling of nearshore processes which affect coastal engineering 
projects. This effort was authorized by Headquarters, U.S. Army Corps of 
Engineers (HQUSACE), under Civil Works Research Work Unit 32484, 
Directionality of Waves in Shallow Water, Coastal Flooding Program. Funds 
were provided through the Coastal Engineering Research Center (CERC), 
U.S. Army Engineer Waterways Experiment Station (WES), under the pro- 
gram management of Ms. Carolyn M. Holmes, CERC. Messrs. John H. 
Lockhart, Jr., John G. Housley, Barry W. Holliday, and David A. Roellig 
were HQUSACE Technical Monitors. 


This report was prepared by Dr. Charles E. Long at CERC’s Field Re- 
search Facility (FRF) in Duck, NC, under the direct supervision of 
Mr. William A. Birkemeier, Chief, FRF, and Mr. Thomas W. Richardson, 
Chief, Engineering Development Division, CERC; and under the general 
supervision of Dr. James R. Houston and Mr. Charles C. Calhoun, Jr., Di- 
rector and Assistant Director, CERC, respectively. 


Data processing software was written by Dr. Joan M. Oltman-Shay of 
Quest Integrated, Inc., Seattle, WA. Physical maintenance of the directional 
gauge was done by the FRF dive team consisting of Messrs. Birkemeier, Mi- 
chael W. Leffler, H. Carl Miller, Eugene W. Bichner, and Brian L. 
Scarborough. Gauge calibration was maintained by Mr. Kent K. Hathaway, 
FRF. Data acquisition and storage were done by Mr. Clifford F. Baron, 
FRF. Considerable assistance in data processing was provided by students 
contracted through the Cooperative Education Program: Ms. Wendy L. Smith 
of Old Dominion University, Mmes. Juliana Admadji and Janna Pemberton of 
Florida Institute of Technology, and Ms. Judy Roughton of College of the 
Albemarle, Elizabeth City, NC. Dr. R. T. Guza, Scripps Institution of 
Oceanography, and Mr. Lockhart provided several useful suggestions after 
reviewing a preliminary draft of this manuscript. The efforts of all of these 
individuals are appreciated. 


At the time of publication of this report, Director of WES was Dr. Robert 
W. Whalin. Commander was COL Bruce K. Howard, EN. 


1 Introduction 


Ocean wind waves are among the dominant forcing mechanisms for near- 
shore dynamic processes. Of the broad suite of possible wave regimes, those 
associated with high energy are most important for engineering design. High- 
energy waves produce the most rapid changes in natural or renourished beach 
morphologies and are most likely, by definition, to exceed the strength or 
elevation of a coastal protective structure, resulting in costly damage or exces- 
sive overtopping. Complete knowledge of the types of extremely energetic 
natural wind wave fields that are incident on the nearshore is thus rather 
critical for efficient design of coastal engineering projects. 


A complete wave field definition describes the way in which wave energy 
is distributed in both frequency and direction. Because of the difficulty in 
making high-resolution measurements of wave energy distributions, early 
engineering guidance treated a given sea state as consisting of a train of waves 
at a single frequency traveling in a single direction. While most wave fields 
are more complicated than this, there persists much guidance based on unidi- 
rectional, monochromatic waves, as evidenced in the Shore Protection Manual 
(SPM 1984). Advances in time series measurement and efficient analysis in 
the last three decades have provided considerable knowledge of the frequency 
distributions of wave energy, enabling the development of good theoretical 
models and improved physical models of wave processes with a consequent 
improvement in engineering guidance. 


Detailed measurements of wave direction are more difficult. Only within 
the last decade have practical methods of directional wave measurement been 
available, and these are still expensive and logistically demanding. Until very 
recently, there were no long-term observations from which to extract the 
directional character of nearshore seas with any statistical confidence. To 
alleviate this lack of knowledge, the Field Research Facility (FRF) of the U.S. 
Army Engineer Waterways Experiment Station, Coastal Engineering Research 
Center, deployed a high-resolution directional wave gauge at its site near 
Duck, NC, in September 1986. Measurements covering a 5-year period, 
lasting through August 1991, enable considerable clarification of the direc- 
tional nature of incident wind waves at this site. 


The purpose of this report is to elucidate, by illustration, the nature of the 


frequency and directional distributions of wind wave energy during storm 
events occurring during the 5-year observation period. These results are 
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intended to assist numerical modelers in developing and testing theories of 
wave transformation across a continental shelf. Results also can be used by 
numerical and physical modelers to establish realistic seaward wind wave 
boundary conditions for models of nearshore processes. Illustrations are con- 
tained primarily in Appendix A. Chapters 2 and 3 of this report describe the 
experiment site, gauges, sampling pattern, storm definition, characterizing 
parameters, and format of displayed results. Common features of storm wave 
frequency-direction spectra are described, and some conjecture about unusual 
spectral properties, along with the appropriate path of future research, is 
presented in Chapter 4. 
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2 Measurement Scheme 


Measurement Site 


The FRF is located on the Outer Banks of North Carolina, just north of the 
village of Duck (Figure 1). The site is representative of many coastal sites 
characterized by a broad continental shelf, relatively mild bottom slope, and 
reasonably uniform longshore topography. The mild sloping bottom substan- 
tially dissipates wind waves in the nearshore so that attention can be focused 
on waves arriving from seaward. The uniform longshore topography allows 
the use of a spatial array of sensors to detect wave direction without contami- 
nation by a longshore inhomogeneity in the incident wave field. The facility, 
described by Birkemeier et al. (1985), has a permanent staff and computing 
facilities that are used to monitor, maintain, and error check the instrumenta- 
tion on a day-to-day basis, and so minimize data losses from the myriad prob- 
lems that can plague field data collection. 


Regional and local topography 


At the large environmental scale, the coastline in the vicinity of the FRF is 
nearly straight for several tens of kilometers north and south. It is oriented 
such that the shore-normal direction (facing seaward) is very nearly 70 deg 
from true north (70 deg T). Waves and onshore winds can approach this site 
along an easterly 180-deg arc from 340 to 160 deg T. The adjacent continen- 
tal shelf is approximately 100 km wide, being narrower to the south where 
Cape Hatteras juts eastward. The direction of nearest approach of the 100-m 
isobath (the shelf break) is 10 deg to 15 deg south of east and is about 80 km 
distant from the FRF. A characteristic bottom slope for the shelf is 1 m/km, 
but, at a finer scale, the bathymetry is marked by numerous features of 1- to 
10-km horizontal scales and 10-m vertical scales scattered irregularly across 
the shelf. Waves with periods near 10 sec begin to shoal at depths of 100 m, 
so that considerable refraction of wind waves propagating over the shelf is 
expected (U.S. Army Engineer District, Wilmington 1980). 


Within a few kilometers of the FRF, the bathymetry is more regular. A 
complex bar system exists within about 300 m of the shore (Birkemeier 1984), 
and waves and currents have created some irregular bathymetry of the 600-m- 
long research pier (Miller, Birkemeier, and DeWall 1983). Away from these 
regions, isobaths are nearly shore-parallel. The bottom slope is nearly uni- 
form at about 0.5 m/km and is reasonably dynamically stable, as indicated by 
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Figure 1. FRF location and offshore bathymetry 


repeated surveys (Howd and Birkemeier 1987). Figure 2 illustrates the 


bathymetry over a region 600 m north to 600 m south of the pier and extend- 
ing about 900 m offshore. 


The site is subject to a variety of climates that give rise to a broad diversi- 
ty in frequency and directionally distributed wave energy. Typical wind and 
wave climate for the period covered by this report can be roughly classified as 
four basic types (Miller et al. 1988; Leffler et al. 1989, 1990, 1991, 1992, 
1993). First are low-wind situations, where the wave field is characteristically 
low-frequency swell propagating in from the open Atlantic Ocean. These 
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Figure 2. FRF nearshore bathymetry and instrument location 


waves typically arrive from the east-to-southeast, and are normally present 
between storm events. 


The second and most common type of wave-generating winds are from 


northeasters associated with atmospheric frontal passages. Fronts typically 
arrive from the west and have a northeast-southwest orientation. Active local 
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wave generation occurs after a frontal passage when strong winds arise in the 
northeast quadrant. A third type of storm is an extratropical depression that 
typically forms south of the FRF along the south Atlantic coast. These storms 
tend to migrate northward along the coast, creating high winds and waves 
from the east to northeast. Less frequent but more intense wave generators 
are hurricanes, the fourth type of disturbance. These originate in the equatori- 
al Atlantic, migrate west to northwest, and can make landfall on the gulf and 
Atlantic coasts with very destructive winds and waves. No hurricanes made 
landfall at the FRF during the course of this study, but several passed by in 
the deeper Atlantic, creating a distinctive local field of high, long-period 
waves narrowly distributed in frequency and direction. 


It is noted that the regional topography and climate limit the generalization 
of results obtained at the FRF to sites having similar characteristics. For 
example, one would not expect nearshore dynamics of sites on the Great 
Lakes, gulf, or west coasts of the United States to behave in the same way be- 
cause their climates and bathymetries are different. Nonetheless, fundamental 
physical processes such as nonlinear wave behavior, wave breaking, and air- 
sea interaction, for which insights are gained by studies at the FRF, will be 
more universally applicable; the limitations are not too restrictive, and gener- 
alizations are justified. 


Instrumentation 


Data from two gauges are presented in this report. One gauge consists of 
a set of near-bottom-mounted pressure gauges, which together constitute a 
high-resolution directional wave gauge. This gauge provides estimates of the 
frequency-direction spectra of wind waves which are the basis of this report. 
The other gauge is an anemometer that gives a measure of local wind speed 
and direction. It provides a characterization of the larger wind field that 
generates some of the wave fields measured with the directional wave gauge. 
These gauges are described here briefly for completeness and because their 
properties affect the way data are interpreted. 


High-resolution directional wave gauge 


The directional wave gauge consists of two fundamental parts: (a) an array 
of gauges that sense sea surface displacement, and (b) an analysis algorithm 
that converts the data to an estimate of S(/,@), the frequency-direction spec- 
trum, where f is cyclic wave frequency and 6 is wave direction.’ The sen- 
sors are a Set of nine near-bottom-mounted pressure gauges, laid out in a 
straight line along the 8-m isobath as shown in Figure 2. Spacing between the 
gauges follows the guidance given by Davis and Regier (1977) for the design 
of wave direction detectors. Pressure time series data from each gauge are 
converted to units of equivalent seawater depth, Fourier transformed, and 
converted to sea surface displacement Fourier transforms using linear wave 


For convenience, symbols and abbreviations are listed in the Notation (Appendix B). 
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theory. For each pair of gauges, raw cross-spectral estimates are then com- 
puted, ensemble averaged, and band averaged to create a smoothed cross- 
spectrum for each gauge pair. For each smoothed frequency, a matrix is 
formed of the cross-spectra for all possible pairs of gauges. This matrix is 
then fed to the directional estimation algorithm. 


Wave energy directional distributions are computed using the Iterative 
Maximum Likelihood Estimator (IMLE), which was derived by Pawka (1982, 
1983), and found to be a substantial improvement on the Maximum Likeli- 
hood Estimator (MLE) derived by Davis and Regier (1977). Specifically, the 
algorithm is Equations 1 - 4 in the paper by Pawka (1983) with convergence 
exponent 6 = 1, convergence coefficient y = 5, and a maximum of 30 itera- 
tions. 


In practice, all gauges of the FRF array are not used for all wave frequen- 
cies. Gauges with large separations are used for long, low-frequency waves. 
Closely spaced gauges are used for shorter waves to avoid loss of coherence 
for these waves at large spatial separations. Details of array sub-configura- 
tions and their associated frequencies are given by Long and Oltman-Shay 
(1991). Tests of the FRF array with artificial data indicate a directional reso- 
lution of 5 deg. This means it can resolve the peak of a single mode direc- 
tional distribution to within 5 deg, or distinguish the peaks of a bimodal distri- 
bution when separated by 15 deg or more. 


Anemometer 


The reference anemometer used in this study was mounted atop the FRF 
laboratory building (Figure 2) at an elevation of 19.5 m above mean sea level. 
It has a nominal speed accuracy of +0.5 m/sec and direction accuracy of 
+3 deg. Speed and direction data are collected on separate channels, and 
results from each channel are simply averaged (not vector averaged). The 
proximity of the anemometer to the FRF building is known to cause low speed 
readings on some occasions in very high winds (Long and Hubertz 1988), but 
it is considered sufficiently reliable for the generally moderate winds encoun- 
tered in this study. Though other anemometers have been deployed for short- 
er periods at the FRF, this one was consistently in use for the duration of this 
study. It follows that this gauge provides the most consistent local represen- 
tation of the larger scale wind field. 


Data Sampling Schemes 


With few exceptions, data used in this report were obtained following the 
standard FRF sampling scheme. The pattern is to collect four contiguous 
records of 34-min duration (4,096 data points at a 2-Hz sampling rate) starting 
at 0100, 0700, 1300, and 1900 Eastern Standard Time (EST). When the 
characteristic wave height H_, exceeds 2 m at a reference wave gauge (usual- 
ly gauge #1 of the directional array), additional collections are made starting 
at 0400, 1000, 1600, and 2200 EST. This form of sampling gives good data 
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coverage of high energy conditions, but has a tendency to miss some of the 
details of initial transient conditions at the onset of a storm. 


The anemometer data presented in this report are based on a subset of the 
overall FRF database. Wind speeds and directions are deduced from single 
34-min records with start times of 0100, 0700, 1300, and 1900 EST. These 
particular records are the most reliable, having been checked by hand for use 
in the FRF annual reports (Leffler et al. 1992, 1993). The result is a very 
regular sampling of local winds, with gaps only for system maintenance or 
repair. 


All possible data were processed from the directional wave gauge, so that 
normally there are samples from all the routine start times listed above. Four 
records of 34 min each (2 hr 16 min or 16,384 data points at 2 Hz) for each 
gauge were used in cross-spectral computation to ensure high confidence in 
the resulting cross-spectral matrix that is transformed in the IMLE algorithm. 
Fifteen half-overlapping segments of 2,048 points were processed for ensem- 
ble averages, and 10 subsequent frequency bands were averaged for a final 
resolution bandwidth Af = 0.00977 Hz and a minimum of 160 degrees of 
freedom (excluding any gain in degrees of freedom from the lapped segments) 
in the final cross-spectral estimates. While this method gives very good re- 
sults in estimates of S(f,@), it has the effect of averaging all wave processes 
for the 2 hr 16 min records so that highly transient processes are smoothed 
somewhat. Results are not invalid, just rather strongly filtered. Shorter 
records can be used in future work if brief transients are important. 


The bulk of S(@) results were obtained by the method just described, but 
there were, on occasion, minor variations in sampling times and rates. Dur- 
ing the SUPERDUCK experiment in September and October 1986 (Crowson 
et al. 1988), data collection start times were adjusted to conform to high and 
low tides, and so did not conform to the FRF normal collection pattern. 
Furthermore, a total of seven 34-min records were obtained in each collection. 
To maximize the number of S(f,@) observations during SUPERDUCK, re- 
cords 1 - 4 from a collection were processed, and then records 4 - 7 were 
processed, so that record 4 was used twice, but observations were still almost 
independent. Time series of parameters from these spectra thus appear rather 
oddly spaced, with pairs of samples appearing with | hr 42 min separations. 
During the DELILAH experiment of October 1990 (Birkemeier et al. 1991), 
data were collected for an additional 34 min and the sampling rate for ail data 
was 4 Hz instead of 2 Hz. Thus, a collection had 10 records of 4,096 data 
points (about 17 min) each. For the data used here, wind speeds and direc- 
tions during DELILAH were averaged for only 17 min instead of the usual 
34 min. Directional array data were still processed in 2 hr 16 min total re- 
cord lengths, but, in smoothing, 20 bands were averaged instead of 10, result- 
ing in the same final resolution bandwidth (but more degrees of freedom). 
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General Working Data Set 


The final working data set from which to extract storm data consists of 
6,762 observations of S(f,@) spanning the period from September 1986 to 
August 1991, or about 5 years. Two major gaps in data collection occurred 
during this time, affecting the seasonal coverage of storms and reducing their 
total number. After SUPERDUCK, full time collection of directional gauge 
data did not resume until late January 1987, so there is a gap in the time line 
of collections from November 1986 to January 1987. In December 1987, a 
fishing boat trawled a net along the line of the linear array of gauges, destroy- 
ing about half of them. Parts replacement and repairs could not be effected 
until August 1988, so there is another gap from December 1987 to August 
1988. Aside from these two major gaps, the record is substantially complete, 
excepting short periods for routine maintenance and calibration. 
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3 Storm Data 


Event Definition 


To extract storm data from the larger working data set, a definition of what 
is meant by a storm is necessary. For the purposes of this report, a storm is 


defined using the characteristic wave height H,, and the length of time it 


exceeds a certain level. To isolate the largest events, it was required that H__, 
exceed 2 m for times at least 24 hr apart, and not drop below | m in the 
interim. Wave conditions are normally relatively quiescent just before a 
storm, generally being low-amplitude, long-period swell from the southeast 
quadrant. To ensure at least one prestorm condition, data for a given event 
were extracted from the main database from the beginning of the day contain- 
ing a prestorm condition to the end of the day during which H_ first drops 
below 1 m. 


Figure 3 shows a time series of H_, from the directional gauge for all col- 
lections in the main database. Storm events that satisfy the above definition 
are identified by letter or number. Their nominal durations are marked by 
horizontal black bars. In all, 29 storm events were found and are referenced 
by the letters A to Z and the numbers | to 3. Note that Event 3 follows 
Event E. 


Table 1 lists the storm events shown in Figure 3, along with the start and 
end dates defined by the criteria described above. Also shown are the maxi- 
mum H_ observed during the storm and a nominal duration, defined as the 
time elapsed from the first time H_,, exceeds 2 m to the last time it drops 
below 2 m. Linear interpolation between measured H_, was used to obtain 
the exceedence times. Based on the definition of a storm described above, all 
durations should exceed 24 hr. Several do not, which illustrates the pitfalls of 
searching a 6,762-point database by eye. The shorter storms are retained here 
because they do not differ from the longer storms except by their durations. 


Climatological Parameters 


Several characterizing parameters are used in this report to describe a 
storm evolution. These parameters are given names and symbols common to 
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Storm Events 
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Figure 3. Time series of H,,, with storm events identified 
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= 
Table 1 
Storm Event Identification 
oe 

H duration 
Event Start End m hr 
A 9 Oct 86 14 Oct 86 3.09 57.01 
B 17 Oct 86 21 Oct 86 2.38 23.67 
(S 15 Feb 87 20 Feb 87 4.03 40.07 
D 9 Mar 87 18 Mar 87 4.22 161.84 
E 22 Mar 87 26 Mar 87 2.68 28.15 
3 11 Oct 87 17 Oct 87 3.11 85.29 
F 10 Nov 87 12 Nov 87 2.45 27.04 
G 26 Nov 87 1 Dec 87 22) 86.13 
H 14 Dec 88 17 Dec 88 2.31 27.34 
| 3 Jan 89 6 Jan 89 2.19 8.68 
J 20 Jan 89 27 Jan 89 3.00 78.70 
K 16 Feb 89 20 Feb 89 3.19 54.43 
L 22 Feb 89 26 Feb 89 4.32 52.32 
M 2 Mar 89 5 Mar 89 2.11 24.45 
N 6 Mar 89 14 Mar 89 4.23 108.62 
ie) 21 Mar 89 | 25 Mar 89 ASS) 30.08 
P 3 Sep 89 11 Sep 89 2.59 138.10 
Q 18 Sep 89 28 Sep 89 2.68 199.53 
R 22 Oct 89 1 Nov 89 2.63 38.02 
Ss 7 Dec 89 11 Dec 89 4.17 61.76 
T 21 Dec 89 26 Dec 89 4.68 91.00 
U 6 Mar 90 9 Mar 90 2.96 22.25 
Vv 25 Oct 90 28 Oct 90 4.70 44.31 
WwW 17 Nov 90 21 Nov 90 2.51 42.70 
xX 7 Dec 90 10 Dec 90 2.41 24.07 
Y 7 Jan 91 14 Jan 91 3.40 111.84 
Z 18 Apr 91 21 Apr 91 3.38 51.00 
1 18 May 91 21 May 91 2.82 33.22 
2 23 Jun 91 28 Jun 91 2.66 21.05 


wave analysis, but because there are several ways to compute each with the 
data at hand, they are defined specifically here. The parameter set consists of 


characteristic wave height H_, characteristic wave period I characteristic 
wave direction 6,, representative directional spread of wave energy 66,, and 
water depth d. 


Characteristic wave height H_ is defined, as usual, as four times the 
standard deviation of sea surface displacement, based on variances contained 
in the frequency spectrum S(f) between the low- and high-frequency cutoffs. 
In this analysis, S(f) is computed as the average of the surface-corrected 
frequency spectra from each of the nine gauges in the linear array. This does 
not increase degrees of freedom (because all gauges are measuring the same 
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wave field), but it does reduce instrument noise. All analysis is done in dis- 

crete frequency space, and results are defined in terms of N = 28 frequency 

bands from f, = 0.054 Hz (the low-frequency bound) to f, = 0.319 Hz (the 
high-frequency bound) in steps of Af = 0.00977 Hz. The discrete frequency 
spectrum is given the symbol S(f,) so that characteristic wave height is 


1-4/3 siya (1) 


Characteristic wave period, or peak period 7, is defined as the inverse of 
the peak frequency f/ deduced from the frequency spectrum; i.e., the discrete 
frequency at which S(f) is maximum. This is the conventional definition of 
peak period, and is used here to be consistent with other wave studies. Thus, 
iS Mie 


There are a number of ways to define peak direction 6,. The one used 
here is the peak of the discrete frequency-direction spectrum S(f,,6,) at the 
peak frequency; i.e., the discrete direction 6, at which S(f,.8,,) is a maxi- 
mum. Discrete directions 6, range from 6, = -90 deg to 6,, = @,, = 90 deg 


in A@ = 2-deg increments. These angles represent the directions from which 
wave energy arrives relative to shore-normal, with -90 deg being waves trav- 
elling parallel to the beach from the south, 0 deg being waves arriving along 
the shore-normal axis, and +90 deg being waves travelling parallel to the 
beach from the north. A constraint on linear array analysis is that it cannot 
detect waves travelling seaward, so all waves are assumed to arrive from the 
seaward semicircle of azimuths. 


A characteristic directional spread of wave energy 60, also is a very useful 
parameter for these data. It is also defined from S(f,,6) as being the angle 
that subtends the central half of the energy in the directional distribution at the 
peak frequency. Symbolically, an integral function /(@) can be defined as 


-(° So) 
I(f,8) = | cs wire (2) 
where /(f,-90°) = 0 and 1(f,90°) = 1 identically because the integral of 
S(f,@) with respect to direction defines S(f). The arc that contains the central 
half of the energy at a given frequency begins at the angle 6,, for which 
I(f,6,.) = 0.25 and ends at the angle @,, for which /(f,6,.) = 0.75. If 6,, 

and 6,, are determined at f = f,, the peak directional spread is simply 


58, = 6,5 - Oss (3) 


This parameter is robust in that it can be computed for any finite directional 
distribution and is meaningful in that it gives a small number for narrowly 


ls 


Chapter 3 Storm Data 


14 


distributed energy, and a large number for broadly distributed energy. Note 
that its use here differs from that of Long and Oltman-Shay (1991). Long and 
Oltman-Shay (1991) looked at the overall spread of S(f,@). Examined here is 
the energy spread only at f = f,, the most energetic frequency spectral line. 


Water depths d reported here, and used for tide stage determination and 
surface correction of near-bottom pressure signals, are the medians of the 
mean depths determined from the nine pressure gauges for the duration of one 
collection. This was necessary because, while the pressure gauges are accu- 
rate at wind wave frequencies, they are susceptible to long-term drifts in their 
mean signals due to environmental changes and biological activity within the 
gauge sensing volumes. Use of median depths as representative depths for all 
the gauges eliminates the largest deviate mean depths. Comparison of result- 
ing mean depths with the FRF tide gauge indicates an accuracy of about 
0.1 m, which is sufficient for use here. 


Data Presentation 


The primary result of this study is the set of figures contained in Appendix 
A. For each of the events listed in Table 1, there are two types of display. 
One page contains time series of climatological parameters associated with the 
event. Following the climatology page for each event, subsequent pages 
contain sequences of frequency-direction spectral plots for select times within 
the event. Page A3 contains an index of events, showing the page in Appen- 
dix A where climatological parameters are displayed. 


Climatological parameters displayed for the duration of each event are 
those described above, namely H_, ie d, 66, 6. wind direction, and wind 
speed. The reader is referred to the parametric plot for Event A on page A4 
for an example. Wind direction and peak wave direction 6, have been trans- 
formed to the same coordinate system (directions from which winds or waves 
arrive in degrees counterclockwise from shore-normal, as discussed earlier), 
and are drawn on the same set of axes. Plot symbols will then be collocated 
if wind and peak wave directions are aligned. 


Because some events contain a large number of observations of S(f,), it is 
somewhat impractical logistically to display all observations for all events. 
Consequently, only certain select cases are illustrated in pages following the 
parametric plots. A prestorm condition, all observations up to and near the 
peak H_ levels in a storm, and several observations during storm decay are 
displayed for all events. Shaded symbols in the graph of H_ for a given 


event indicate which observations of S(f,@) are plotted. For some storms of 
short duration, all observations are plotted. 


Sequences of observed frequency-direction spectra, selected as noted, are 
illustrated on pages following the parametric plots. For example, selected 
cases from Event A, indicated by the shaded symbols on page A4, are shown 
as sequences of plots on pages AS and A6. For each observation there are 
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four entities displayed: (a) a three-dimensional plot of the surface represent- 
ing S(f,6), (b) a linear graph of the frequency spectrum S(f), (c) a linear 
graph of the direction spectrum S(6) (defined below), and (d) a contour plot 
of S(f,@). The same vertical scale is used for all three-dimensional plots of 
S(f,6) within an event, so the relative importance of each observation can be 
assessed. Prestorm and storm decay conditions thus show low relief, and 
peak conditions tend to have the largest relief. 


The linear plot of S(f) is shown in the right panel of the three-dimensional 
box, with its abscissa parallel to the frequency axis of the S(,@) plot. The 
direction spectrum S(@) is the directional analog of the frequency spectrum 
and shows the directional distribution of total energy. Computationally, it is 
found from 


S(6) = ( S(f,6) df (4) 


and has units of square meters per degree. The analog with S(f) can be seen 
from the definition of the frequency spectrum as 


Si) = | Ki S(f,) dé (5) 


with units of square meters per Hertz. The left panel of the three-dimensional 
box contains a linear graph of 5(@). Both S(f) and S(@) are drawn to constant 
scales derived from the maxima of these entities within an event. 


The contour plot of 5(f{6) is shown below the three-dimensional box for 
each observation. Contour intervals are tenths of the maximum value of 
S(f,6) of all observations in an event. Consequently, spectra having relatively 
low energy levels have both low relief in the three-dimensional plots and few, 
if any, lines in corresponding contour plots. Observations with higher energy 
levels have correspondingly more contour lines. 


Because all four plotted entities within a given event are drawn to constant 
scales within an event, axis annotation is rather redundant. It is also space 
consuming. Consequently, axis annotation is drawn only for the upper left- 
hand set of plots on each page illustrating sequences of S(f,@). Axis annota- 
tion for all other plot axes within an event is identical. 
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4 Discussion 


The primary purpose of this report is to elucidate the directional distri- 
butions of nearshore wind wave energy in the context of evolving storms. 
With the illustrations in Appendix A for guidance, some fundamental proper- 
ties of storm wave fields can be identified. 


A primary question is whether a storm wave field can, in any sense, be 
reasonably approximated as unidirectional and monochromatic, or, perhaps, as 
corresponding to simple models of deep-ocean wave fields. A quick perusal 
of Appendix A shows, in a few observations, that a unidirectional, monochro- 
matic model may be appropriate for some applications (examine, e.g., Event 
R and conditions near the peaks of Events V and W), but such conditions are 
relatively rare. Furthermore, narrow spreads of naturally occurring spectral 
energy in frequency and direction do not mean the waves are unidirectional 
and monochromatic. Real ocean waves are invariably groupy and subject to 
nonlinear interactions. These dynamically important properties can not be 
represented by unidirectional, monochromatic wave models, and the utility of 
such models is thereby somewhat restricted. 


The wave fields of Appendix A do not appear to conform to models and 
observations of deep-ocean wave fields, either. For example, papers by 
Mitsuyasu et al. (1975) and Mitsuyasu and Uji (1989) indicate that S(4@) in 
the deep ocean is characterized by narrower directional spreads near the spec- 
tral peak frequency than in the spectral tails. Spectra shown in Appendix A 
appear generally to have maximum spreads near spectral peak frequencies. 
These observations mean that neither a single-wave-train model nor deep- 
ocean wave behavior is appropriate to characterize nearshore storm waves. 


While the spectra of storm wave fields at the FRF are complicated, they 
are not without a remarkable similarity in structure that can be described, at 
least qualitatively. In this part of the report, these structural similarities are 
discussed, and some of the consequent research questions are posited. 


Of the 29 storm events identified in this report, 25 have S(@) with com- 
mon structural features (described qualitatively below) through the progress of 
a storm. Four events, E, P, Q, and R, do not follow this pattern. Of these 
four, waves in Events P and Q are associated with Hurricanes Gabrielle and 
Hugo, respectively, both of which passed several hundred kilometers away 


Chapter 4 Discussion 


from the FRF. Well-dispersed, low-frequency wave trains with very narrow 
directional spreads are associated with these events and are characteristic of 
waves from distant sources. Events E and R may be natural deviates from the 
common pattern of the remaining 25 events. 


Storm waves driven by wind systems identified in the FRF annual reports 
(Miller et al. 1988, Leffler et al. 1989, 1990, 1991, 1992, 1993) as Canadian 
high-pressure systems, frontal passages, or low-pressure cells forming west to 
south of the FRF have frequency-direction spectra with many common fea- 
tures. Most of these features are illustrated in Event N (pages A47 to ASO), 
which exemplifies the set of events. The progression of these features as the 
storm evolves is as follows: 


a. There is usually a prestorm wave field consisting of low energy, long 
period swell with moderate spread arriving from shore-normal to slight- 
ly south of shore-normal. For Event N, this is the first case illustrated 
on the top row of page A48. 


b. During the primary growth stage, at the initial onset of winds, high- 
frequency, locally generated waves appear. Peak period usually drops 
from the pre-storm condition. Peak direction shifts nearly to align with 
wind, but not exactly (bottom-mounted pressure gauges can detect only 
intermediate to shallow water waves, so some refraction towards shore- 
normal will take place even in short fetches). Peak directional spreads 
are of order 10 deg to 20 deg, the smallest of any time in the storm. 
For Event N, this is the second case illustrated on the top row of page 
A48. 


c. In the secondary growth stage (remainder of the top row and first two 
cases of the second row of plots on page A48 for Event N), peak period 
increases, and the frequency spectrum tends to conform to one of the 
conventional nearshore models, such as TMA (Bouws et al. 1985) or 
FRF (Miller and Vincent 1990). It is the directional distribution within 
S(f,@) that is most curious. As energy at the peak period increases, so 
does the characteristic spread, perhaps indicating continued local growth 
plus additional energy radiating in from growth further offshore. More- 
over, an additional lobe of energy appears at the peak frequency, 
aligned more nearly with shore-normal or, for Event N, well south of 
shore-normal. Furthermore, a secondary peak of high-frequency energy 
(identified as a ridge in three-dimensional plots of S(/,@) or a second 
tail off the main lobe of energy in contour plots) appears, making a 
strongly bimodal high-frequency tail to the spectrum. The second 
energy ridge is aligned nearer to shore-normal, and ranges from 20 deg 
to 40 deg offset from the initial growth stage spectrum. 


d. In the advanced growth stage (last two cases on the second row and first 
two cases on the bottom row of plots on page A48 for Event N), the 
spectral peak continues to move to lower frequencies and becomes 
centered at very nearly shore-normal. The spread is rather large, being 
of order 40 deg to 50 deg, and the secondary lobe of low-frequency 
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energy appears to be absorbed by the primary peak. The structure of 
the bimodal, high-frequency tail remains, and its energy increases. 


e. As storm waves reach their peak wave energy, the directionally broad, 
energetic, low-frequency spectral peak and directionally bimodal, high- 
frequency tail of the spectra remain. If driving winds persist, this 
structure can continue relatively unchanged for some time. In Event N, 
for example, wind speed remains nearly constant at 16 to 18 m/sec, and 
wind direction persists at 60 deg to 80 deg from shore-normal for more 
than 2 days (see page A47). The sequence of observations on pages 
A48 to A50, showing S(f,6) from 1300 on 7 Mar 89 to 1900 on 9 Mar 
89, illustrate a structure that changes very little, suggesting that some 
equilibrium condition has been attained. 


f. In the waning stages of a storm (the final plots on page A50 for spectra 
of Event N), winds usually decrease in speed and change direction. 
Active local wave generation abates, and the high-frequency part of the 
spectrum decays. Low-frequency waves that were generated offshore 
continue to arrive, still with a broad directional spread. 


g. In the end stage, conditions generally return to the prestorm spectrum 
of low-frequency, low-energy swell, having broad directional spread 
and being centered near the shore-normal direction. 


The above description was derived specifically from Event N, but the other 
24 storm events have the same features to a greater or lesser degree. Even 
the two dissimilar storms, Events E and R mentioned above, have some of the 
same features. The ubiquitousness of these spectral patterns indicates a high 
degree of consistency in storm wave evolution at the FRF. The unexpected 
features of these spectra, low-frequency side lobe of energy in the secondary 
growth stage, as well as a directionally broad, low-frequency peak coupled 
with a bimodal high-frequency tail in the most energetic stages, require some 
explanation so that guidance appropriate to engineering projects can be de- 
rived. 


There are a number of lines of approach for explaining these observations. 
Several of them require work that is beyond the scope of this report, but the 
nature of the necessary investigation can be mentioned here. A primitive 


question is whether the observed S(f,@) are real or some artifact of the inter- 
action of the true wave field and the IMLE algorithm as used with the FRF 
linear array of gauges. If the wave energy is all propagating onshore, there is 
no known reason to suspect the spectral estimation process. Tests using syn- 
thetic data indicate the system is well capable of resolving the broad direction- 
al spreads and 20-deg-separated bimodal distributions observed here. If there 
is wave energy reflected from nearshore, the directional estimation system 
used in this report cannot resolve it, and the reflected energy will appear to be 
coming from offshore. This question can be addressed with a more recently 
installed 360-deg array, but that analysis has not been completed at this writ- 
ing. Nonetheless, it is unlikely that reflections will account for the present 
observations because the spectral pattern recurs at many tidal stages and over 
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many seasons, wherein nearshore reflectivity would be expected to vary con- 
siderably. 


Another possibility is that linear wave theory is not adequate during high- 
energy conditions when the wave field may become saturated with highly 
nonlinear, actively breaking waves. While this may be true for the extreme 
storm of Event T, it is unlikely to be true for the very marginally energetic 
storm of Event I, which was inadvertently included in this analysis, and had 
an H_, that exceeded 2 m for less than 9 hr. The common spectral pattern 
described above also occurs in Event I (see pages A33 and A34), so it is 
improbable that a saturated wave field accounts for these observations. 


Two other possible causes for these results are a common structure in the 
overall driving wind fields or the net effect of the broad, shallow bathymetry 
of the continental shelf adjacent to the FRF. The anemometer used in this 
report is only a crude indicator of the larger wind systems that generate waves 
observed at the FRF. A possible line of future research to examine the larger 
systems would be to apply a wind field model like that used 
in the CERC Wave Information Studies (Resio, Vincent, and Corson 1982) to 
the storm events examined in this report. 


In conjunction with this possible study would be the inclusion of a reason- 
ably accurate bathymetry for the United States mid-Atlantic continental shelf. 
The shelf break, at depths of 60 to 100 m, is 80 to 100 km distant from the 
FRF (Figure 1). Waves with periods of 8 to Il sec (common periods for 
waves observed at the FRF) begin to be influenced by the bottom at these 
depths. Consequently, considerable wave refraction is expected for intermedi- 
ate and long waves generated over the shelf and in the deeper Atlantic by 
nonuniform wind fields associated with atmospheric fronts and low pressure 
systems. This condition may account for the directionally broad, low-frequen- 
cy peaks and side lobes in observed storm spectra. However, a quantitative 
assessment of these effects should be made, either with a numerical model or 
with more extensive observations. 


It is possible that the directionally bimodal, high-frequency spectra may be 
due to a hummock or linear bar just offshore of the FRF array that causes 
strong refraction for high-frequency waves. This hypothesis can be tested 
with a detailed bathymetric survey that extends beyond the area routinely mea- 
sured by the FRF staff (generally, the area shown in Figure 2). Such a sur- 
vey has not yet been done, but the results presented herein should provide 
strong motivation for such work. There is weak evidence against a local 
bathymetric effect in the present data set. In the later stages of Events A, Y, 
and 2, winds shifted to the southeast with just enough force to generate waves 
from that quadrant. Though rather faintly shown in the figures (see pages A6, 
A87, and A96), bimodal spectral structures appear in the southeast quadrant in 
forms roughly symmetric to those observed in the northeast quadrant. It is 
somewhat unlikely that a local bathymetric feature would induce this kind of 


symmetry in S(f,6). 
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Another possible explanation for the present observations is nonlinear 
interactions among wave trains during their generation and propagation to the 
FRF directional gauge. Various interaction mechanisms have been proposed 
by the WAMDI Group (1988), Freilich, Guza, and Elgar (1990), Resio and 
Perrie (1991), and others. Testing these hypotheses requires more detailed 
spatial measurements of S(f,6) than are currently available. However, an 
experiment can and should be designed for this purpose if the current results 
are to be fully explained. 


By way of guidance for model testing and engineering design, a representa- 
tive example of S(f,@) for late stage growth and peak storm conditions is 
shown in Figure 4. It is simply an average of the 19 observations from the 
relatively stationary peak storm conditions of Event N. It contains the prima- 
ry features noted above and represents a nearly stationary dynamic process 
that continued for 56 hr. It should be noted that although the high-frequency 
tail has a spectral density less than the spectral peak, the total energy in the 
part of the tail that is bimodal (found from the area under the frequency spec- 
trum over the high-frequency range of frequencies) is roughly 25 percent to 
50 percent of the total energy. Because the bimodal tail occurs at peak energy 
conditions, there is considerable energy in the tail region. 


The key point is that a wave generation and transformation model at conti- 
nental shelf and larger scales clearly must have a result that looks qualitatively 
like Figure 4 to represent conditions at the FRF with reasonable verity. Addi- 
tionally, nearshore process models of the FRF will be improved if a spectrum 
like that in Figure 4 is used to describe the wave boundary condition outside 
the surf zone in high-energy cases. 
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Figure 4. A characteristic storm spectrum from the FRF 
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5 Conclusion 


This report elucidates the nature of frequency-direction spectra evolving 
during storm events occurring at the FRF. The primary result is a series of 
graphic presentations of spectra measured with a high-resolution directional 
wave gauge in an 8-m depth approximately 900 m offshore. Data from a 
5-year observation period have been searched for storm events defined >v a 
characteristic wave height in excess of 2 m at two points at least 24 hr “art 
in the same storm system. A total of 29 events have been isolated for closer 
scrutiny. Time series plots of climatic parameters and sequences of plots of 


S(f,@) are used to reveal the evolutionary nature of storm wave fields. 


The data indicate that there are two basic types of events: (a) high waves 
from strong but distant sources, and (b) localized wave generation from Cana- 
dian high-pressure systems, frontal passages, and atmospheric low-pressure 
cells passing the FRF from the west or south. Of the 29 events identified, 27 
are of this latter type. Of the 27 localized events, 25 demonstrated very com- 
mon patterns of spectral evolution. The patterns include initial localized 
waves aligned nearly with the local wind, directional broadening of spectral 
peaks as the wave fields grow, and a bifurcation of high-frequency spectral 
tails into two modes: (a) one nearly aligned with the wind, and (b) another, 
20 deg to 40 deg away, aligned nearer to the shore-normal direction. One- 
quarter to one-half of the total energy is contained in the bimodal tails, a 
consideration of consequence because such spectra occur at the peaks of 
storms, and because the beating of such spectral waves against the shoreface 
can affect such nearshore processes as the patterns of beach erosion and infra- 
gravity wave generation. 


The main conclusions here are that storm waves occur at the FRF with 
curious spectral shapes, and that accounting for the processes that give rise to 
those shapes is important. It is suggested that a more detailed coverage of 
wind, wave, and bathymetric observations over 100-km scales from the FRF 
is warranted. It is certain that nearshore processes driven by waves having 
such elaborately structured spectra will be different from estimates generated 
by over-simplified wave models (typically unidirectional and monochromatic) 
as even now still haunts engineering guidance. There is good reason to be- 
lieve the spectra reported here are real and not some artifact of the analysis 
method used. Use of these results for testing wave generation and transforma- 
tion models at scales of the continental shelf width, and as seaward wave 
boundary conditions for nearshore process models, will certainly yield more 
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realistic, meaningful, and accurate representations of nature, and thereby an 
improved coastal engineering science. 
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Frequency-Direction Spectral Plots: Event A 
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Climatological Parameters 
Event B: 17 Oct 86 to 21 Oct 86 
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Frequency-Direction Spectral Plots: Event B 
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Climatological Parameters 
Event C: 15 Feb 87 to 20 Feb 87 
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Frequency-Direction Spectral Plots: Event C 
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Climatological Parameters 
Event D: 9 Mar 87 to 18 Mar 87 
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Frequency-Direction Spectral Plots: Event D 
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Frequency-Direction Spectral Plots: Event D 
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Event E: 22 Mar 87 to 26 Mar 87 
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Frequency-Direction Spectral Plots: Event E 
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Frequency-Direction Spectral Plots: Event E 


0400, 24 Mar 87 
Hmo = 2.01m 0700, 24 Mar 87 1000, 24 Mar 87 1300, 24 Mar 87 
Hmo = 1.96 m Hmo = 2.08m Hmo = 1.88 m 


BAQ\ KKNDOO et Ide) 


90 Cy) 
@ (deg) 


1900, 24 Mar 87 0100, 25 Mar 87 0700, 25 Mar 87 1900, 25 Mar 87 
Hmo = 1.71m Hmo = 1.34 m Hmo = 1.17 m Hmo = 1.05 m 


Thess 


0100, 26 Mar 87 
Hmo = 1.10 m 


0700, 26 Mar 87 
Hmo = 0.99 m 


1300, 26 Mar 87 
Hmo = 1.04 m 


1900, 26 Mar 87 
Hmo = 1.10 m 


A19 


Appendix A Event Time Series and Spectral Sequences 


Climatological Parameters 
Events 3: ll OcteS7ato u7Octy Si 


sopOFoFOn, 


oo0000o Qa 


eo Foegeeee 
ogoo00o000000 


depth (m) 


direction (deg) 


wind speed (m/sec) 


A20 


Appendix A Event Time Series and Spectral Sequences 


Frequency-Direction Spectral Plots: Event 3 
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Event F: 10 Nov 87 to 12 Nov 87 
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Frequency-Direction Spectral Plots: Event F 
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Frequency-Direction Spectral Plots: Event F 
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Event G: 26 Nov 87 to 1 Dec 87 
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Frequency-Direction Spectral Plots: Event G 
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Frequency-Direction Spectral Plots: Event H 
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Frequency-Direction Spectral Plots: Event | 
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Event J: 20 Jan 89 to 27 Jan 89 
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Frequency-Direction Spectral Plots: Event J 
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Frequency-Direction Spectral Plots: Event K 
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Frequency-Direction Spectral Plots: Event N 
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Frequency-Direction Spectral Plots: Event O 
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Frequency-Direction Spectral Plots: Event Q 
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Frequency-Direction Spectral Plots: Event Q 
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Frequency-Direction Spectral Plots: Event R 
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0100, 7 Dec 90 
Hmo = 0.42 m 0700, 7 Dec 90 1300, 7 Dec 90 1900, 7 Dec 90 
Hmo = 0.41m Hmo = 0.42 m Hmo = 0.44 m 


iy 


BAO) N00 wt [ded 


CYS) 


90 t) 
6 (deg) 


0100, 8 Dec 90 0700, 8 Dec 90 1300, 8 Dec 90 1600, 8 Dec 90 
Hmo = 1.08 m Hmo = 1.71m Hmo = 2.24 m Hmo = 2.41 m 


1900, 8 Dec 90 2200, 8 Dec 90 0100, 9 Dec 90 0400, 9 Dec 90 
Hmo = 2.14 m Hmo @ 2.16 m Hmo = 2.17 m Hmo = 2.34 m 


A83 


Appendix A Event Time Series and Spectral Sequences 


Frequency-Direction Spectral Plots: Event X 
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Frequency-Direction Spectral Plots: Event Y 
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Frequency-Direction Spectral Plots: Event Z 
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Frequency-Direction Spectral Plots: Event Z 
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Climatological Parameters 
Event 1: 18 May 91 to 21 May 91 
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Frequency-Direction Spectral Plots: Event 1 
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Frequency-Direction Spectral Plots: Event 1 
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Frequency-Direction Spectral Plots: Event 2 
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Frequency-Direction Spectral Plots: Event 2 
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Appendix B Notation 


Water depth 

Cyclic frequency 

Cyclic frequency of n” component wave 
Spectral peak frequency 

High-frequency boundary for a discrete spectrum 
Low-frequency boundary for a discrete spectrum 
Spectrum-based characteristic wave height 
Direction integral of directional distribution function 
Iterative Maximum Likelihood Estimation 

Index integer 

Upper limit of indices denoted m 

Maximum Likelihood Estimation 

Index integer 

Upper limit of indices denoted n 

Spectral density 

Frequency spectral density 


Frequency spectral density at n” discrete frequency 
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S86) Frequency-direction spectral density 


S(f,,6,,) Frequency-direction spectral density at n” discrete 
frequency and m” discrete direction 


S(6) Direction spectral density 

i Characteristic wave period derived from spectral peak 
6 Convergence exponent in IMLE method 
Y Convergence coefficient in IMLE method 

66, Directional spread of wave energy at spectral peak 

frequency 

Af Finite frequency increment 
6 Wave direction 

@). Element m of a discrete set of wave directions 

o., Upper bound of discrete direction coordinate 

6, Peak wave direction 

6, Lower bound of discrete direction coordinate 

Oh Angle subtending 25 percent of wave energy 

6, Angle subtending 75 percent of wave energy 
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